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Abstract Solar wind protons observed by the MESSENGER spacecraft in orbit about Mercury exhibit 
signatures of precipitation loss to Mercury's surface. We apply proton-reflection magnetometry to sense 
Mercury's surface magnetic field intensity in the planet's northern and southern hemispheres. The results are 
consistent with a dipole field offset to the north and show that the technique may be used to resolve 
regional-scale fields at the surface. The proton loss cones indicate persistent ion precipitation to the 
surface in the northern magnetospheric cusp region and in the southern hemisphere at low nightside 
latitudes. The latter observation implies that most of the surface in Mercury's southern hemisphere is 
continuously bombarded by plasma, in contrast with the premise that the global magnetic field largely 
protects the planetary surface from the solar wind. 


1. Introduction 

A remarkable feature of Mercury's weak, internal magnetic field, indicated by orbital observations, is a 
-480 km northward offset of the magnetic equator from the geographic equator [Anderson et al., 201 1, 2012; 
Johnson et al., 201 2]. The low magnetic field strength and the northward offset provide constraints on 
Mercury's enigmatic dynamo mechanism [Stanley and Clatzmaier, 2010] and lead to direct interactions 
between the solar wind and the planet's surface. The weak magnetic field allows precipitation of solar 
wind plasma to the surface in the region of the northern magnetospheric cusp. Because the surface magnetic 
field strength at any southern latitude is predicted to be weaker than that at the corresponding northern 
latitude, enhanced ion sputtering and space weathering in the southern hemisphere are possible but to date 
have been expected to occur mainly on the dayside, at high latitudes in the southern cusp region. In this study 
we apply proton-reflection magnetometry, adapted from electron reflectometry [Lin etal., 1 976; Mitchell et al., 
2001; Lillis et al., 2008], to determine Mercury's surface magnetic field strength in both hemispheres and 
measure particle precipitation to the surface. 

Electron reflectometry (ER) has been used at the Moon [Lin et al., 1 998], Mars [Lillis et al., 2008], and Ganymede 
[Williams et al., 1997] to sense remotely the magnetic field strength at the surface. ER depends on the 
magnetic mirroring effect, that is, the reflection of electrons by the gradient in the magnetic field strength 
along convergent field lines. Electrons that would mirror below the surface are lost, and the flux of reflected 
electrons exhibits a sharp drop at the pitch angle a (the angle between the particle velocity and the local 
magnetic field direction) that corresponds to mirroring at the surface. The in situ magnetic field together with 
the pitch angle of the last reflected electrons, the cutoff pitch angle, indicates the surface magnetic field 
strength. This technique has not yet been applied using protons. 

At Mercury, protons with energies of 0.3-10 keV are regularly detected inside the magnetosphere [Zurbuchen 
et al., 201 1] by the Fast Imaging Plasma Spectrometer (FIPS) [Andrews et al., 2007] on the MErcury Surface, Space 
ENvironment, GEochemistry, and Ranging (MESSENGER) spacecraft. The Energetic Particle Spectrometer [Andrews 
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Figure 1 . (a) Proton counts (green dots) versus zenith angle, 0, and azimuth angle, <(>, in the FIPS reference frame, together 
with pitch angle contours for an example integration period (magenta contour is for a = 90°, green contour is for a = 150°, 
and the contour interval is 10°). The spacecraft longitude, latitude, and altitude are also given, (b) Time-normalization 
weighting factor as a function of a. (c) Derived PAD at this integration period. 

et al., 2007] on MESSENGER detects electrons with energies only above 35 keV and observes substantial fluxes 
too infrequently to use ER. We therefore use FIPS observations of protons within Mercury's magnetosphere 
[ Winslow et al., 201 3], together with the magnetic field observations from MESSENGER'S Magnetometer (MAG) 
[Anderson et al., 2007], to estimate the magnetic field strength at Mercury's surface. 

2. Deriving Pitch Angie Distributions 

For our analyses, we used one Earth year of observations, from 7 June 201 1 to 7 June 201 2. Combining 
observations from MAG and FIPS allowed the calculation of proton pitch angle and the derivation of pitch 
angle distributions (PADs) within the magnetosphere. 

FIPS measures energy per charge (E/q), time of flight, and arrival incidence angle for ions, and it completes 
one scan over the full range of E/q values every 8 s [ Raines et al., 201 1], The sensor has a conical field of 
view (FOV) of approximately 1 An sr, with two symmetric cutouts of 1 5° aperture, one near the 
instrument's symmetry axis and another near the plane perpendicular to the symmetry axis. FIPS pulse 
height analysis (PH A) data specific to proton times of flight are used to create the proton PADs at each 
integration period of FIPS. The raw, three-dimensional PFIA events (Figure 1 a) were divided by the visible 
solid angle per detector pixel and the pixel-dependent microchannel plate (MCP) detector efficiency 
[Raines et al., 2014]. Because only angular structure is of interest for these accumulations, proton events 
from all E/q steps were added together to improve the signal-to-noise ratio. 

After rotating the vector magnetic field measured by MAG into the FIPS reference frame, pitch angles for 
all proton counts were computed. The solid-angle-weighted and MCP-normalized proton counts that fall 
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b) 


within 10° pitch angle contours were then 
summed into 10°-wide pitch angle bins, centered 
on pitch angles from 5° to 175°. 

The accumulated proton events were divided by a 
time-normalization factor that takes into account 
the FIPS visibility of each magnetic pitch angle bin 
[Raines et at., 2014]. This normalization is a function 
of the FIPS orientation with respect to the magnetic 
field direction and FOV obstructions at each 
integration time (Figure 1 b). We thus arrive at a PAD 
at every integration time, for which the proton 
counts have been properly weighted in order to 
take into account all the viewing limitations 
(Figure 1c). Errors assigned to each pitch angle bin 
center incorporate counting statistics and the 
fraction of proton gyrophase angles that was visible 
in the FOV for each pitch angle bin. 

Because of the limited field of view of FIPS, the 
full proton PAD from a = 0° to 180° is not visible at 
any one integration time. To build PADs spanning 
et = 0° to 180° at times when the spacecraft was 
over the same location in Mercury solar orbital 
(MSO) coordinates (for which +X is sunward, +Z is 
northward, and -t-Tcompletes the right-handed 
system), we combined and averaged PADs from 
times when the incident population was 
observed with PADs obtained when the reflected 
population was visible. We focused on regions 
where the highest proton counts were detected 
by FIPS. This analysis thus yields Mercury's 
average long-wavelength (e.g., dipole) surface 
magnetic field structure but does not resolve 
shorter-wavelength structure. 

To connect the incident and reflected 
populations from different integration periods, 
we normalized each individual PAD by the 
weighted proton count observed at a = 85° or 95°, 
depending on which of the two pitch angles was 
visible in the distribution at that time, to approximate the average incident population. This methodology 
follows that applied to electron reflectometry measurements at the Moon [Halekas et at., 2001]. We then 
computed the weighted averages of the individual PADs, where the inverse of the error at each pitch angle bin 
center on the PADs was assigned as the weight. The error assigned to the final averaged distribution is 1 
standard error in all the proton counts averaged at each pitch angle bin center. 



a (degrees) 

Figure 2. (a) PAD for Mercury's northern cusp, (b) PAD for the low- 
latitude southern hemisphere nightside. Average proton counts 
are in red with standard errors. The black curve shows the diffusion 
model fit to the reflected portion of the distribution; the fit 
uncertainty is in gray. The yellow shading indicates the loss cone; 
the black error bar shows the uncertainty in a c . 


3. Surface Field Estimates 

Averaged PADs were derived in the northern magnetospheric cusp region (geographic latitude A > 60°N), 
as well as in regions of high proton flux on the nightside at low latitudes (0° <A < 30°S) in the southern 
hemisphere. Although high proton fluxes to the dayside, in the southern cusp region, are expected, 
MESSENGER is unable to measure these, because its eccentric orbit allows the spacecraft to be inside the 
magnetosphere on the dayside only over northern latitudes. The averaged PAD for the northern cusp, 
which includes 485 reflected and 185 incident-population scans, is shown in Figure 2a. The incident 
population is approximately independent of pitch angle (i.e., isotropic), consistent with protons entering 
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along open field lines from the 
magnetosheath. However, there is a 
loss cone in the fluxes of reflected 
protons, from a~1 20° to 1 80°, evidence 
for incident protons having been "lost" 
to Mercury's surface. Relative to the 
sharp cutoff pitch angle observed in ER 
at the Moon [Mitchell et at., 2008], the 
edge of the proton loss cone is 
smoothed over -30° in a. We attribute 
this more gradual transition to the 
combined effects of the FIPS angular 
resolution of ~1 5° and pitch angle 
diffusion from wave-particle scattering 
in the cusp. Broadband magnetic field 
fluctuations of 0.001 -10 Hz frequency 
are consistently observed in this region 
[Anderson et at., 2013] and would pitch 
angle scatter protons. 

We estimated the loss cone angle 
consistent with diffusive scattering by 
fitting solutions to the diffusion equation 
to the loss cones of our pitch angle 
distributions. We solved the one- 
dimensional diffusion equation 


8u(a, t) d 2 u(a, f) 

' = LJa ~ 


8t 


8a 2 


(D 
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Figure 3. (a) Absolute values of misfit of the diffusion model to the PAD 
in the northern cusp, as a function of a c and D a t. The bold contour marks 
the upper bound on the misfit level from which the errors on the fit 
parameters were obtained, (b) Same as in Figure 3a but for the southern 
hemisphere PAD. The ratio of the bold contour misfit level to the mini- 
mum misfit is ~7 in Figure 3a and ~3 in Figure 3b. 

the average maximum and minimum weighted proton counts in 
given by 


where u is the proton count, f is the time 
over which diffusion occurs in the PAD, 
and D a is the diffusion coefficient, with a 
step function initial condition: 


{ Ci, a < a c 

( 2 ) 

c 2 , a > a c 


where c ^ and c 2 are constants set by 
the PAD. The boundary conditions were 


8u( 0, t) 8 u(7t, t) 


8a 


da 


= 0 . 


(3) 


The solution to equation (1) is 


u(a , t) = ^2 B n cos (na)e n2D ° r , 


(4) 


where 


B n = 


J Ci cos(na)da+ J c 2 cos (na)dc 


(5) 


We fit equations (4) and (5) to our loss cones allowing the cutoff pitch angle, a c , and D a t to vary. We used a 
grid search method that minimized the median absolute deviation between the model and the observations. 
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Figure 4. (a) Stereographic projection (looking at the planet from above the north pole) of surface magnetic foot-point 
locations of reflected proton observations in the northern cusp. The Sun is to the right, (b) Corresponding plot for the 
southern hemisphere (looking at the planet from above the south pole). 

Figure 3 shows the contour plots of the absolute value of this misfit measure for the models best fit to the 
averaged PADs. We established the upper and lower bounds on the parameters a c and D„f by identifying an 
allowable upper bound on the misfit, corresponding to a 95% confidence limit (bold contours in Figure 3). 
The upper and lower bounds on a c were identified as the locations of the intersection of a horizontal cut 
(passing through the minimum misfit) with the bold contour. A corresponding vertical cut yielded the limits 
on the D a t parameter. The bounds on the best fit model, shown by the gray shaded regions of Figure 2, were 
determined from the diffusion curves corresponding to the upper and lower limits for a c and D a t. 

From the fit to the northern cusp PAD, we obtain a cutoff pitch angle of 121°±3°, which together with the 
measured average magnetic field strength at spacecraft altitudes (<550 km), B 0 = 302.4 ± 53.0 nT, implies a surface 
field strength of B s = 41 2 ± 98 nT from the relation B s = 6 0 /sin 2 (a c ). The uncertainty accounts for the standard error 
in the fit value of a c and also the standard deviation of B 0 , computed from all the intervals in the average PAD. 

We mapped the observation locations in the cusp to the surface by tracing field lines to the surface given 
Mercury's offset dipole field. The observation altitudes ranged from 282 km to 549 km, with a mean of 414 km. 
Cusp observations traced to the surface extended 15.6° in latitude and 7.5 h in local time and were centered 
on noon at 76.4°N latitude (Figure 4a). 

We also find high proton fluxes in the latitudinal band 0° <2< 30°S on the nightside, with a clear loss cone 
signature in the PAD (Figure 2b), although with larger uncertainties than for the northern cusp (Figure 2a). 
In the southern hemisphere, observations as far south as possible are desirable for observing the long- 
wavelength structure of the magnetic field. Flowever, because of MESSENGER'S eccentric orbit, we are 
restricted to observations north of ~30°S latitude. In this averaged PAD, we included 128 scans in the 
reflected population and 315 scans in the incident population. The similar errors on most of the proton fluxes 
in Figure 2b, despite the higher number of observations included for the incident side, result from large 
standard deviations in the incoming fluxes. The diffusion model that best fits the loss cone gives a c = 43°+ 2 3 , 
and this value together with B 0 = 52.5 ± 1 4.8 nT at the spacecraft altitude corresponds to a surface field of 
B s = 1 1 3+|, nT. The observation altitudes were higher than those in the northern cusp, ranging between 
1 1 60 and 1 980 km, with a mean of 1 535 km. The mapped surface locations span 23°S < X < 34°S, with a mean 
of 27.8°S, and local times spanning the nightside from 16 h to 5.3 h, centered on 23.5 h (Figure 4b). The 
apparent secondary loss cone in the incident population in Figure 2b implies that these observations 
correspond to closed field lines on the nightside. 

4. Discussion and Conclusions 

The estimates of surface magnetic field strength are compared with predictions from the best fit time- 
averaged magnetospheric model of Johnson et al. [2012] in Figure 5 and Table 1. The results from proton- 
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Figure 5. Surface magnetic field strength from the best fit magnetospheric model in MSO coordinates, Mollweide projection. 
The offset of the magnetic equator from the geographic equator and the magnetopause and magnetotail fields cause 
departures from the field of a centered dipole alone. The heavy black lines mark the approximate outline of regions 
sampled by proton-reflection magnetometry. 



s 


reflection magnetometry are markedly lower than the model magnetic field. Such a difference is expected, 
however, because the model is for a vacuum magnetic field, whereas our PADs demonstrate the presence of 
plasma that extends to the surface. The plasma generates a diamagnetic field, which reduces the surface field 
below the vacuum model prediction [ Korth et al., 201 1, 2012; Winslow et a!., 2012]. 

Fortunately, the proton data provide the information required to estimate the diamagnetic effect. The 
particle flux at the surface in the northern cusp can be determined from the loss cone size and mean proton 
temperature, 7" p , and density, n p , in the cusp, yielding the proton pressure P p = n p kT p , where k is Boltzmann's 
constant. The typical proton density (/i p = 30 cirT 3 ) and temperature (T p ~ 12 MK) derived from FIPS 
observations in the cusp [ Raines et al., 2014] are for an isotropic particle distribution; the anisotropy 
associated with the loss cone produces an underestimate of plasma density [Gershman et al., 2013]. Taking 
this anisotropy into account yields a surface flux of 3.7 x 10 12 particles nrT 2 s -1 , which is approximately in 
agreement with the flux determined by Winslow et al., [201 2] for the cusp region. For our southern 
hemisphere average, we have n p = 5 citT 3 and T p ~ 20 MK, which give a flux of 4.4 x 10 11 particles mr 2 s~\ 
The diamagnetic effect reflects particle motions only in the direction perpendicular to the magnetic field, so 
we used the perpendicular particle pressure to calculate the diamagnetic field at the surface. We find an 
average diamagnetic effect, AB P | asma , in the northern cusp of 86 ± 1 1 nT and AB P | asma = 44 ± 6 nT in the 
southern hemisphere region. Uncertainties in these values represent the limits on the diamagnetic field 
derived from the pressure at the spacecraft altitude and that at the surface. After accounting for this effect, 
our magnetic field estimates agree with the model predictions in the southern hemisphere and are within 8% 
(our upper bound) of the model in the northern cusp (Table 1). 

The validity of the offset dipole can be tested by estimating the ratio between the surface magnetic field 
strength in the northern cusp and that in the southern hemisphere. From the mean observation locations, we 
find that the ratio from proton-reflection magnetometry, corrected for the diamagnetic effect, is 3.2+ 3 °, in 
agreement with the result of 3.3 for an offset dipole field alone. The full magnetospheric model [Johnson 
et al., 201 2] yields a surface field ratio of 2.6 to 3.6 for the latitudinal and longitudinal extent of our 


Table 1. Surface Magnetic Field Strength From Proton-Reflection Magnetometry Compared With Magnetospheric 


Model Predictions 
Latitude 3 (MSO Surface) 

Local Time 3 (MSO) 

B P r (nT) b 

^Plasma ( n T) 

^PR ^Plasma ( n "0 

^Model ( n T) 

76.4°N 

27.8°S 

1 2 h 
23.5 h 

412 ±98 

1 1 3 +87 
1 lo -61 

—86 ± 1 1 
-44 ±6 

498 ± 99 
157^ 

647 

224 


a Central latitudes and local times of the two PADs. 

“Surface magnetic field estimates from proton-reflection magnetometry. 

“Magnetospheric model predictions for comparison with estimates from proton reflectometry corrected for 
plasma effects. 
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observation locations, with a mean value of 2.9. A centered dipole field alone gives a ratio of 1.5 for the 
latitudes sampled and is not in agreement with our inferred surface field ratio. 

Proton-reflection magnetometry thus provides independent confirmation of the offset dipole nature of Mercury's 
magnetic field, although weaker terms of higher degree and order in a multipole expansion of the field are not 
ruled out by these results. Significantly, proton-reflection magnetometry demonstrates that persistent proton 
precipitation to the surface occurs on the nightside of Mercury's low-latitude southern hemisphere. This result 
implies that such precipitation may also be occurring at mid-southern latitudes on the nightside, where closed 
field lines reach the surface, because the northward offset of the magnetic equator results in weak surface field 
strengths (not more than 25% above the value from proton-reflection magnetometry) everywhere in the 
southern hemisphere. Together with the proximity of the magnetopause to the surface and the large cusp region 
on the dayside southern hemisphere [Anderson et at., 2011, 2012; Winslow et a!., 2012], this finding suggests 
that most of Mercury's southern hemisphere surface may be continuously bombarded by plasma. Such 
precipitation in the southern hemisphere implies that space weathering is not confined to the cusp and may thus 
show limited latitudinal variation. Because of MESSENGER'S eccentric orbit, proton-reflection magnetometry on 
the nightside at higher southern latitudes is not feasible, but observations from the ESA-JAXA BepiColombo 
mission [Benkhoffet at., 2010], with a less eccentric orbit that will provide low-altitude observations in both 
hemispheres, may provide an opportunity to probe the surface field more extensively. 
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